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ABSTRACT 

We report the discovery of a massive (MpSim= 13.02 ± 0.64 Mj; total mass 13.25 ± 0.64 Mj), 
large (1.95 ± 0.16 Rj) planet in a transiting, eccentric orbit (e = 0.260 ± 0.017) around a lO"^ 
magnitude F5V star in the constellation Camelopardalis. We designate the planet X0-3b, and the 
star XO-3, also known as GSC 03727-01064. The orbital period of X0-3b is 3. 1915426 ± 0.00014 
days. XO-3 lacks a trigonometric distance; we estimate its distance to be 260±23 pc. The radius 
of XO-3 is 2.13±0.21 Rq, its mass is 1.41±0.08 Mq, its t;sim = 18.54 ± 0.17 km s~\ and its 
metallicity is [Fe/H] = —0.177 ± 0.027. This system is unusual for a number of reasons. X0-3b 
is one of the most massive planets discovered around any star for which the orbital period is 
less than 10 days. The mass is near the deuterium burning limit of 13 Mj, which is a proposed 
boundary between planets and brown dwarfs. Although Burrows et al. (2001) propose that 
formation in a disk or formation in the interstellar medium in a manner similar to stars is a 
more logical way to differentiate planets and brown dwarfs, our current observations arc not 
adequate to address this distinction. X0-3b is also unusual in that its eccentricity is large given 
its relatively short orbital period. Both the planetary radius and the inclination are functions 
of the spectroscopically determined stellar radius. Analysis of the transit light curve of X0-3b 
suggests that the spectroscopically derived parameters may be over estimated. Though relatively 
noisy, the light curves favor a smaller radius in order to better match the steepness of the ingress 
and egress. The light curve fits imply a planetary radius of 1.25±0.15 Rj, which would correspond 
to a mass of 12.03 ± 0.46 Mj. A precise trigonometric parallax measurement or a very accurate 
light curve is needed to resolve the uncertainty in the planetary mass and radius. 

Subject headings: binaries: eclipsing - planetary systems - stars: individual (GSC 03727-01064) - tech- 
niques: photometric - techniques: radial velocities 
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1. Introduction 

There are now over 200 extrasolar planets 
known (http://exoplanets.eu/), and most of these 
have been discovered using the radial velocity 
technique. As a result, most of these systems 
yield the minimum mass of the planet (MpSim), 
its orbital semi-major axis, and properties of the 
central star. The data to date have dramatically 
changed our appreciation of the diversity of plane- 
tary systems that can form, and have particularly 
focussed the community's attention on the role 
that planetary migration plays in the planet for- 
mation process (see review by Papaloizou et al. 
2007). As the precision and duration of radial 
velocity surveys increases, lower mass planets and 
planets in longer orbits continue to be found. As 
a result, it has been suggested that the fraction of 
stars hosting planets may be as high as 50% (see 
Udry et al. 2007). Therefore, it appears that the 
planet formation process is relatively efhcient. 

The generally accepted model of giant planet 
formation is that of core nucleated accretion (e.g.. 
Pollack et al. 1996; Bodenheimer et al. 2000; 
Hubickyj et al. 2005) in which gas giant plan- 
ets first form a several 10) earth mass core 
of solids via non-elastic collisions in a disk fol- 
lowed by the runaway accretion of gas from the 
disk once the mass of the core is sufficient to ex- 
ert a strong gravitational influence. Such a model 
naturally predicts a positive correlation between 
the metallicity of the disk (and parent star) and 
the ease with which planets can form: there are 
more building blocks to form the solid core in 
higher metallicity cases. The expected correlation 
has been suspected for some time (e.g., Gonza- 
lez 1997; Santos et al. 2003), and has now been 
demonstrated in an unbiased way (e.g., Santos et 
al. 2004; Fischer & Valenti 2005) and has been 
taken as evidence in support of the core accretion 
model of planet formation. The competing model, 
gravitational instability in a massive circumstellar 
disk (e.g.. Boss 1997, 2000; Mayer et al. 2002) 
does not predict such a relationship. Livio and 
Pringle (2003) showed that lower disk metallicity 



03016, Guarcino (FR) Italy 

^^Bellatrix Observatory, Via Madonna dc Loco, 47 03023 
Ceccano (FR) Italy 

i*CBA Belgium Observatory, Walhostraat lA, B-3401 
Landen, Belgium 



reduces the efficiency of Type II migration which 
may account for some of the increased probability 
for finding a planet around more metal rich stars. 
There is some observational evidence suggesting 
a higher frequency for planets at smaller separa- 
tions around metal-rich stars which may support 
such a metallicity-migration relationship (Jones 
2004; Sozzetti 2004). As a result, the planet- 
metallicity correlation may (at least in part) be 
due to a metallicity-migration relationship and 
may not simply signify a metallicity-formation re- 
lationship. 

Additional observations of extrasolar planets 
are necessary to distinguish between the compet- 
ing modes of planet formation and the conditions 
which affect planet migration scenarios. For ex- 
ample, the core accretion model requires several 
Myr to form Jovian planets in a disk (Inaba et 
al. 2003), while gravitational instabilities could 
potentially form these planets around very young 
stars. Jovian planets are suspected around a few 
^ 1 Myr young stars (e.g., CoKu Tau/4 ~ Forrest 
et al. 2004; GQ Lup - Neuhauser et al. 2005); 
however, clear detection and firm mass determina- 
tions from techniques such as radial velocity vari- 
ations remain elusive on such young stars. Tran- 
siting extrasolar planets offer the opportunity to 
determine both the mass and radius of the planet, 
and hence the planet's density. Such planets al- 
low us to better constrain the variety of extraso- 
lar planet properties and offer another potential 
way to distinguish between planet formation sce- 
narios. For example, the high core mass derived 
for the transiting Saturnian mass planet in orbit 
around HD 149026 (Sato et al. 2005) has been 
interpreted as strong support for the core accre- 
tion model. There are now ~ 20 known transiting 
extrasolar planets (http: / /exoplanets.eu/) . Sev- 
eral of these planets have unexpectedly large radii 
(see Figure 3 of Bakos et al. 2007 and discussion 
therein), perhaps suggesting additional heating of 
these planets other than that from irradiation by 
the star they orbit. Additional effort is needed to 
understand the structure of extrasolar planets and 
their origins. 

Here, we report the discovery of the third tran- 
siting planet from the XO Project. The XO 
Project aims to find planets transiting stars suffi- 
ciently bright to enable interesting follow up stud- 
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ies (McCuUough et al. 2005). Q This planet 
presents some very interesting properties: it has a 
mass of Mpsim= 13.02 ± 0.64 Mj and is in a short 
period (3.2 d), eccentric orbit (e = 0.260 ± 0.017) 
around the apparently single F5V star GSC 03727- 
01064. X0-3b has a radius of 1.95±0.16 Rj, mak- 
ing it substantially larger in both radius and mass 
than most of the other reported transiting plan- 
ets. In §2 we discuss the observations leading to 
the discovery of X0-3b. In §3 we present our anal- 
ysis of these observations to determine the stellar 
and planetary properties. In §4 we give a discus- 
sion of X0-3b in the context of other extrasolar 
planets and current planet formation models, and 
§5 summarizes our conclusions. 

2. Observations 

2.1. XO Project Photometry 

McCuUough et al. (2005) describe the instru- 
mentation, operation, analysis, and preliminary 
results of the XO Project. In summary, the XO 
observatory monitored tens of thousands of bright 
{V < 12) stars twice every ten minutes on clear 
nights for more than 2 months per season of vis- 
ibility for each particular star, over the period 
September 2003 to September 2005. From our 
analysis of more than 3000 observations per star, 
we identified XO-3 (Figure [T]) as one of dozens 
of stars with light curves suggestive of a transit- 
ing planet. With the XO cameras on Haleakala, 
we observed four transits or segments of tran- 
sits of XO-3 in 2003 and three in 2004, on Ju- 
lian dates 2452934, 2452937, 2452998, 2453001, 
2453301, 2453320, and 2453352. Table H] provides 
a sample of the photometry for XO-3 from the 
XO cameras. The full table is available in the on- 
line edition. From the survey photometry of XO-3 
(Figure [2|) , which has a nominal standard devia- 
tion of 0.8% or 8 mmag per observation, we de- 
termined a preliminary light curve and ephemeris, 
which we used to schedule observations of higher 
quality with other telescopes, as described in the 
next subsections. 



This paper includes data taken on the Haleakala summit 
maintained by the University of Hawaii, the Lowell Ob- 
servatory, the Hobby-Eberly Telescope, the McDonald Ob- 
servatory of The University of Texas at Austin, and five 
backyard observatories. 



2.2. Additional Photometry 

As outlined by McCuUough and Burke (2006), 
once an interesting candidate is detedcted in the 
XO photometry, the candidate is released to an 
Extended Team (E.T.) of sophisticated amateur 
astronomers for additional observation. In 2006 
September through November and 2007 January 
through March, transit events of XO-3 were ob- 
served by members of the E.T. from a total of 4 
backyard observatories. Table [T] also provides E.T. 
photometry for XO-3. For the E.T. light curves, 
the median differential magnitude out of transit 
provides the flux normalization and the standard 
deviation out of transit provides the uncertainty 
in the measurements. A fifth E.T. observatory 
was used to obtain all-sky photometry for XO-3. 
These observatories are equipped with telescopes 
of aperture 0.3 m. The telescopes, equipped 
with CCD detectors, are suitable for obtaining 
light curves sufficient to confirm the nature of the 
transit and obtain good timing information. A 
network of such telescopes is well suited to observe 
candidates with known positions and ephemerides. 
Here, we use two such telescopes in North Amer- 
ica and two in Europe. Unlike the case of XO-lb 
(McCuUough et al. 2006) and X0-2b (Burke et al. 
2007) , we were not able to schedule observing time 
on a larger 1-m) telescope to obtain additional 
photometry - the E.T. photometry is the best we 
have obtained to date on X0-3b. 

We estimate all-sky photometric B, V, Rq, and 
Iq magnitudes for XO-3 and several nearby refer- 
ence stars (Figure [T] and Table [2]) calibrated using 
a total of 6 Landolt areas (Landolt 1992). Using 
the Hereford Arizona Observatory (E.T. member 
BLG) 0.36-meter telescope on photometric nights 
2006 October 27 and 2006 November 5, we mea- 
sured the fluxes of 48 and 18 Landolt stars, re- 
spectively at an air mass similar to that for XO-3 
and established the zero points of the instrumental 
magnitudes and transformation equations for the 
color corrections for each fllter and the CCD. The 
derived magnitudes for XO-3 differed on the two 
dates by 0.03 magnitudes or less in all 4 colors. 
The B,V,Rc, and Ic absolute photometric accu- 
racies are 0.03, 0.03, 0.03, and 0.05 mag r.m.s., 
including both the formal error and an estimated 
systematic error. The Tycho magnitudes for XO-3 
listed in Table [3] transform (via Table 2 of Bessel 
2000) to Johnson V = 9.86, i.e. 0.06 mag (2-cr) 
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fainter than our estimate. 

2.3. Spectroscopy 

In order to measure the orbital elements of the 
system, and in particular the mass ratio, as well as 
to determine the characteristics of the host star, 
we obtained spectra of XO-3 with two-dimensional 
cross-dispersed echelle spectrographs (TuU et al. 
1995; Tull 1998) at the coude focus of the 2.7-m 
Harlan J. Smith (HJS) telescope and via a fiber 
optic cable on the 11-m Hobby-Eberly Telescope 
(HET). Both telescopes are located at McDonald 
Observatory. The HJS spectra were obtained in 
a traditionally scheduled manner, while the HET 
spectra were obtained in queue scheduled mode. 
An iodine gas cell is used on the HET to provide 
the wavelength reference for velocity determina- 
tion. While an iodine cell is available for HJS tele- 
scope, it is not a facility instrument and the obser- 
vations of XO-3 collected here were done as a part 
of a radial velocity survey of young stars (Huerta 
2007; Huerta et al. 2007) for which very high (few 
m s~^) velocity precision is not required. Wave- 
length and resulting velocity calibration is accom- 
plished by taking thorium-argon reference lamp 
spectra before and after each observation of XO-3 
at the HJS telescope. At the HJS we obtained one 
spectrum per night, and at the HET we obtained 
or two spectra per night. At both telescopes, 
the spectral resolution was R = A/AA ~ 60000, 
and data were obtained on a total of 21 nights. 
The two-dimensional echelle spectra were reduced 
using IDL procedures described in Hinkle et al. 
(2000) which include bias subtraction, flat field- 
ing using a quartz lamp spectrum, and optimal 
extraction of the data. Table [4] gives a log of the 
spectral observations. 

3. Analysis 
3.1. Ephemeris 

We follow the same procedure used by Burke 
et al. (2007) for X0-2b to refine the ephemeris of 
X0-3b. We adopt the transit midpoint calculated 
for the event of 2006 October 16 as our ephemeris 
zeropoint: 2454025. 3967±0.0021. AU measured 
transit midpoints used to refine the ephemeris are 
reported in Table El To determine the orbital 
period, we minimize the difference between 
the observed transit times and a constant-period 



ephemeris model. Based upon 3 transit events ob- 
served either in multiple passbands by the same 
observer, or observed by different observers, we es- 
timate the I — a uncertainty in the time of the cen- 
ter of an individual transit to be 5 minutes, which 
we adopt as the uncertainty for each of the transits 
observed in a single bandpass by only one observer. 
The best-fit period is 3.1915426±0.00014days. 

3.2. Radial Velocity Measurements 

Using spectra obtained at the HET, we mea- 
sured XO-3's radial velocities with respect to the 
topocentric frame using iodine absorption lines su- 
perposed on the spectra of XO-3. We modeled the 
extracted spectra using a template stellar spec- 
trum and the absorption spectrum of the HET io- 
dine gas cell (Cochran 2000). For the analysis of 
XO-1 (McCuUough et al. 2006) and XO-2 (Burke 
et al. 2007), a high resolution spectrum of the 
Sun and the Earth's atmosphere (Wallace et al. 
1998) was used for the template stellar spectrum. 
We followed this same procedure initially for XO- 
3; however, the resulting radial velocity uncertain- 
ties were substantially higher than we achieved for 
XO-1 and XO-2. We suspect that the higher vsmi 
of XO-3 compared to these other two stars con- 
tributes to this increased uncertainty, but we were 
also concerned that the higher temperature of XO- 
3 relative to these stars (and the Sun) resulted in 
spectral differences large enough to increase the 
uncertainty further. Therefore, we repeated the 
HET radial velocity determinations using a high 
resolution {X/6X ^ 60000) spectrum of the F5V 
star HD 30652 from the SPOCS sample of stars 
(Valenti & Fischer 2005, hereafter VF05) as the 
template stellar spectrum. The effective temper- 
ature of HD 30652 is 6424 K (VF05), within 5 
K of our derived T^jf for XO-3 (see below), and 
z;sim = 16.8 km s^^ which is close to the value of 
18.54 we find for XO-3 below. The spectrum of 
HD 30652 is very similar in appearance to that of 
XO-3. 

Using an implementation of Neldcr and 

Mead's (1965) downhill simplex minimization 
algorithm, "Amoeba," we adjusted parameters of 
our model spectrum to fit the observations. The 
model includes convolution of our model spectra 



^IDL is a software product of ITT Visual Information Solu- 
tions. 
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with a best-fitting Voigt profile to approximate 
the (slightly non-Gaussian) line-spread function 
of the instrument. The free parameters of our 
model are a continuum normalization factor, the 
radial velocity of the star, the radial velocity of 
the iodine lines (which represent instrumental de- 
viations from their expected zero velocity with re- 
spect to the observatory), and an exponent (op- 
tical depth scale factor) that scales the depths of 
the lines as an arbitrary method of adjusting the 
spectrum of HD 30652 to even more closely match 
that of XO-3. Due to the iodine absorption, we 
could not estimate the continuum level by inter- 
polating between local maxima in the spectrum, so 
instead we solved for the continuum iteratively, as 
required to improve the fit between our model and 
the observations. In the manner described above, 
for each ^^15 A section of each individual spectrum 
within the region of the recorded spectrum with 
significant iodine absorption, 5100 — 5700 A, we 
estimated the radial velocity of the star. From the 
approximately normal distribution of the resulting 
radial velocity estimates for each epoch, we calcu- 
lated the mean radial velocity and its uncertainty. 
The l-fj internal errors of the radial velocity mea- 
surements from the HET spectra range from ~ 120 
to ~ 200 m s~^ per epoch. These uncertainties are 
an order of magnitude larger than those achieved 
on XO-1 (~ 15 m s'^) and XO-2 (- 20 m s'^). 
We transformed our measured radial velocities to 
the barycentric frame of the solar system and sub- 
tracted the mean radial velocity and report these 
values in Table HI These velocity measurements 
are shown in Figure U) 

The HJS radial velocity measurement technique 
used here is described in detail in Huerta (2007) 
and Huerta et al. (2007). We summarize it here. 
The radial velocity shift of each observed spec- 
trum is determined by a cross correlation analysis 
of the observed spectrum with respect to a ref- 
erence spectrum. To avoid complications which 
might result from a spectral type mismatch, we use 
one particular HJS epoch (2454137.8215) of XO-3 
as the reference spectrum. Again, wavelength cal- 
ibration is done by averaging the wavelength so- 
lution from Thorium- Argon reference lamps taken 
before and after each stellar observation. A to- 
tal of 15 spectral orders which contain numerous, 
strong stellar lines and no detectable telluric ab- 
sorption lines are used in the cross correlation 



analysis. The radial velocity shift from the 15 
orders are averaged to get the final radial veloc- 
ity shift for each observation, and the standard 
deviation of this mean is computed and adopted 
as the uncertainty of the cross correlation anal- 
ysis. Barycentric radial velocity corrections are 
then applied to the observations to determine ac- 
curate relative velocities for all the HJS spectra of 
XO-3. 

In addition to the uncertainty found above from 
the order-to-order scatter in the HJS radial veloc- 
ity measurements, there is an uncertainty associ- 
ated with the fact that our wavelength calibration 
is based on spectra (Thorium- Argon) that are not 
observed simultaneously with the stellar spectrum 
as is the case for iodine cell observations. In or- 
der to evaluate this uncertainty, a total of 6 stars 
from the Lick Planet Search sample (Nidever et 
al. 2002; Butler et al. 1996; Gumming et al. 
1999) which are known to be stable at the 5 - 
20 m s""'^ level are observed each night and ana- 
lyzed in exactly the same way as we treat the spec- 
tra of XO-3. The internal uncertainty in the ra- 
dial velocity shift measurements for these standard 
stars based on the measured order-to-order scat- 
ter is typically < 20 m s~^. For the 2007 Febru- 
ary observing run (observations 3-10 at the HJS), 
the time series of radial velocity measurements for 
these standard stars show an average standard de- 
viation of 123 m s~^ , with all but one of these stars 
showing a measured standard deviation less than 
126 m s^^. Therefore, we adopt the average value 
of 123 m s^^ as the uncertainty due to the non- 
simultaneous nature of the HJS wavelength cali- 
bration for the 2007 February HJS observations 
of XO-3 and add it in quadrature to the uncer- 
tainty from the order-to-order scatter for each ob- 
servation. For the 2006 October observations, the 
radial velocity standard star data gives an uncer- 
tainty of 132 m s~^ due to the non-simultaneous 
nature of the target and wavelength calibration 
spectra. In addition, Huerta (2007) and Huerta et 
al. (2007) show that there are systematic offsets of 
the order of 100 m s~^ in the radial velocity stan- 
dard star measurements from one observing run to 
the next, and speculate that this is caused by re- 
placing the spectrometer slit plug between observ- 
ing runs (the slit plug is not moved during the ob- 
serving runs). Using the radial velocity standard 
star data, we estimate that there is a systematic 
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shift of 140 ± 33 m s'^ between the 2006 October 
and 2007 February observations. We correct our 
2006 October radial velocity measurements of XO- 
3 by this amount and add this associated uncer- 
tainty in quadrature with the other radial velocity 
uncertainties described above. The measured ra- 
dial velocity shifts and the total uncertainty are 
reported in Table [4] and shown in Figure |4] phased 
to the ephemeris known from the transits. The 
uncertainty for the radial velocity (= m s~^) of 
the reference epoch (2454137.8215) is set equal to 
the lowest uncertainty determined for the other 
HJS relative radial velocity measurements. 

An eccentricity approximately equal to zero is 
expected theoretically for hot Jupiters in ~ 3 day 
orbits (Bodenheimer et al. 2001) and was our ex- 
pectation for X0-3b. Therefore, we expected si- 
nusoidal radial velocity variations for XO-3 with a 
phasing consistent with the transit observations. 
The first two observations obtained with the HJS, 
taken at phase ~ 0.64 and ~ 0.95 are inconsis- 
tent with this expectation, forcing us to consider 
eccentric orbits. Many additional data points are 
required to adequately constrain an eccentric or- 
bit, so we performed intense observing of XO-3 
with both the HJS and the HET in late 2006 and 
early 2007. When fitting the orbit, in addition 
to the radial velocity data, we also use the time 
of mid transit as a constraint in the fitting. Ad- 
ditionally, since the HJS and HET data are on 
different relative velocity scales, we treat as a free 
parameter the offset between these two scales. In 
the orbit fitting, we keep as fixed the orbital period 
and phases determined from the transit ephemeris, 
and treat as free parameters the center of mass 
velocity of the system, the eccentricity, the veloc- 
ity amplitude K for XO-3, the longitude of peri- 
astron, the phase of periastron passage, and the 
offset between the HJS and HET radial velocity 
measurements. We use the nonlinear least squares 
technique of Marquardt (see Bevington & Robin- 
son 1992) to find the best fit parameters for the 
orbit, which are given in Table [6l Figure [4] shows 
the radial velocity data along with the orbital so- 
lution. In this plot, the HET data are shifted (by 
adding 1192.52 m s~^ to the HET radial velocities) 
to match the HJS spectra and the center of mass 
velocity (1002.68 m s~^) of all points is subtracted 
for display purposes, since our center of mass ve- 
locity is based largely on only relative data. (As a 



result of these velocity shifts, the HJS velocities in 
Table m have had 1002.68 m s~^ subtracted from 
them for display in Figure HI and the HET veloc- 
ities in Tabled] have had 189.84 m s~^ added to 
them for display in Figured) Uncertainties in the 
orbital fit parameters are derived by Monte Carlo 
simulation of the data: for 1000 simulations we 
construct fake radial velocity data using the or- 
bital fit and applying Gaussian random noise at a 
level commensurate with the observed data, and 
then fit this model data using the same procedure 
outlined above. Because the orbital fit produces a 
minimum reduced = 0.53, we believe our radial 
velocity uncertainty estimates may be a little too 
large. This can be seen in Table[3]where we report 
the observed minus calculated {O — C) velocities 
of the fit. These values are typically less than our 
derived uncertainties for the radial velocity mea- 
surements. The standard deviation of the HJS 
O — C values is 122 m s~^ compared to a mean 
derived uncertainty of 156 m s~^. For the HET 
data, the standard deviation of the O — C values 
is 90 m s~^compared to mean derived uncertainty 
of 169 m s^^(though this is dominated by a sin- 
gle observation) , indicating that the derived uncer- 
tainties for both data sets are too large. Reducing 
these uncertainties would reduce the uncertainties 
in the fit parameters, but we leave them as given, 
since we can not justify any specific reduction in 
our measured radial velocity uncertainties. From 
the fit, we find the velocity semi- amplitude of the 
orbit of XO-3 is K = 1471±48 m s'^. 

As described below, the vs'mi of XO-3 is rela- 
tively large (18.5 ± 0.2 km s"^), which is substan- 
tially larger than most stars around which planets 
have been discovered (HAT-P-2b is another ex- 
ample of a planet orbiting a star which has a rel- 
atively large vs'mi = 19.8 ± 1.6 km s'-^, Bakos et 
al. 2007). The experience of Mandushev et al. 
(2005) and the similarity between the primary in 
that system with XO-3 gives concern that XO-3 
could be a hierachical triple system in which two 
lower mass stars comprise a short period eclipsing 
binary star which are themselves orbiting a more 
massive, significantly brighter primary in a much 
longer period orbit. In the case of Mandushev 
et al. (2005), the short period binary stars dis- 
play tens of km s~^ velocity shifts, but the strong, 
broad lines of the system almost hide this signal. 
The apparent velocity shift of the system primary 
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is then the result of weaker hnes moving around 
the (stationary) hnes of the primary, producing 
hne profile distortions which are misinterpreted as 
radial velocity variability of the primary. 

Line bisector analysis is a standard technique 
to look for evidence that line distortions are pro- 
ducing false radial velocity variations (e.g., Hatzes 
et al. 1997; Martinez Fiorenzano et al. 2005), and 
have been used successfully to identify false ra- 
dial velocity variations caused by star spots (e.g., 
Queloz et al. 2001; Bouvier et al. 2007; Huerta 
et al. 2007) and hierarchical triple systems (e.g., 
Santos et al. 2002; Mandushev et al. 2005). To 
improve signal-to- noise in the bisector analysis, it 
is common practice to compute the bisector of the 
spectrum cross correlation function (Queloz et al. 
2001; Santos et al. 2002; Mandushev et al. 2005). 
Here, we calculate the line bisector of the total 
cross correlation function computed from all 15 
orders in the HJS spectra used in the determina- 
tion of the radial velocities. We then compute the 
bisector span, which is the difference of the line bi- 
sector at two reference levels. Here we define the 
bisector span, 5b, as the value of the bisector at an 
absolute cross correlation level of 0.15 minus the 
bisector at an absolute cross correlation level of 
0.75. These span values are given in Table ID Be- 
cause the spectral regions of the HET data which 
contain strong stellar lines also contain numerous 
iodine absorption lines, we restrict ourselves to the 
HJS data for the analysis of line bisectors. Previ- 
ous studies which used bisector analysis to reveal 
a false positive in planet searches find a clear cor- 
relation between the measured radial velocity and 
the bisector span (Queloz et al. 2001; Santos et al. 
2002; Mandushev et al. 2005; Bouvier et al. 2007; 
Huerta et al. 2007). In these cases, whether due 
to spots (e.g., Queloz et al. 2001; Bouvier et al. 
2007; Huerta et al. 2007) or due to a hirearchical 
triple system (Santos et al. 2002; Mandushev et 
al. 2005), the full range measured in Sb is approx- 
imately the same as the full range of the measured 
radial velocities. 

On the other hand, Torres et al. (2004) use the 
lack of a correlation between Sb and the radial ve- 
locity to help confirm a planet transiting OGLE- 
TR-56b. In addition, for OGLE-TR-56b the fuU 
range in Sb is less than one third the full range 
in the radial velocities, which Torres et al. (2004) 
argue is additional evidence for the planetary na- 



ture of this object. Figure [5] shows the bisector 
span, Sb, plotted versus the measured radial ve- 
locity for the HJS spectra of XO-3. Also given in 
the figure is the value of the linear correlation co- 
efficient and its associated false alarm probability 
(Bevington & Robinson 1992) for these data. In 
addition to there being no significant correlation 
between S b and the radial velocity, the full range 
of Sb is very small compared to the full range in 
the radial velocity. Figure and Table [3] show that 
the full range in Sb in XO-3 is less than one tenth 
the full range in the radial velocity variations for 
XO-3. 

As a last check that the photometric dimmings 
and radial velocity variations of XO-3 are due to 
a planetary mass companion and not the result of 
a hierarchical triple system with later type stars 
present, we looked for an infrared (IR) excess in 
XO-3. O'Donovan et al. (2006) used a redder 
than expected V — K measurement as one piece 
of evidence to reject the planetary hypothesis for 
the F star GSC 03885-00829. Using the photom- 
etry reported in Table [3l we estimate the color 

V - = 1.01 ± 0.06 for XO-3. Based on the 
effective temperature and gravity found below for 
XO-3, we assign a spectral type of F5V based 
on the calibrations presented in Cox (2000). In 
the color system of Bessell and Brett (1988), an 
F5V star such as XO-3 is expected to have an 
intrinsic {V — Ks)o = 1.10. For comparison to the 
2MASS colors reported in Table[31 we use the color 
transformation relations found in the Explanatory 
Supplement to the 2MASS All Sky Data Release 
and Extended Mission Products by Cutri et al. 

f http://www.ipac. caltech.edu/2mass/releases/allsky/doc/ 1, 
and find that we must add 0.039 mag to the Bessel 
& Brett value in order to estimate {V — Ks)o = 
1.14 on the 2MASS system. Thus, the observed 

V ~ = 1.01 color of XO-3 is 0.13 mag bluer 
(^ 2-cr) than its expected color. The slight bluing 
may be the result of the relatively low metallicity 
we derive (see below), but there certainly is no 
evidence for a near IR excess suggesting a hier- 
archical triple system. The bisector analysis and 
the V — Ks color of XO-3 give us confidence that 
both the photometric dimming and radial velocity 
variations we measure in this star are indeed due 
to a planetary mass companion. 
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3.3. Spectroscopically-Derived Stellar Prop- 
erties and Planetary Mass 

We used the software package SME (Valenti & 
Piskunov 1996) to fit each of the 10 spectra of XO- 
3 from the HJS telescope with synthetic spectra. 
We used the methodology of VF05, including their 
minor corrections to match the Sun and remove 
abundance trends with temperature (negligible in 
this case). Because of gaps between echelle or- 
ders, the HJS spectra are missing the wavelength 
intervals 6000-6123 A, which was included in the 
Valenti & Fischer analysis. These wavelength in- 
tervals are also missing from our extracted HEX 
spectra because the relevant echelle orders span 
the gap between the two detectors. 

We averaged our SME results for the 10 HJS 
spectra, obtaining the parameter values in Table 
[71 Each value in the last column of the table, 
labeled "Precision" because systematic uncertain- 
ties are not included, is the standard deviation 
of the 10 measurements divided by \/9 to yield 
the formal uncertainty in the mean. The median 
value of each derived parameter (not given) dif- 
fers from the mean by less than the uncertainty in 
the mean. The final row in the table gives [Si/Fe], 
which VF05 used as a proxy for alpha-element en- 
richment, when interpolating isochrones. Figure[S] 
shows XO-3's spectrum in the region of the Mg I B 
triplet, which is the dominant spectroscopic con- 
straint on gravity. These three Mg lines also have 
a significant impact on the global [M/H] parame- 
ter, which is used to scale solar abundances for all 
elements other than Na, Si, Ti, Fe, and Ni. 

Following the methodology of Fischer & Valenti 
(2005), we interpolated Yonsei-Yale (Y^) isochrones 
(Demarque et al. 2004) to determine probabil- 
ity distribution functions for the mass, radius, 
gravity, and age of XO-3. The trigonometric 
parallax of XO-3 is unknown, so we initially as- 
sumed distances of 240, 260, and 280 pc with an 
adopted uncertainty of 10 pc in each case (which 
affects the width of the resulting distribution func- 
tions). In order to perform the isochrone analysis, 
we need a V magnitude corrected for extinction 
for XO-3. We measured V = 9.80 ± 0.03 and 
B = 10.25 ± 0.03 for XO-3 (Table E]). These ob- 
servations give a measured B — V = 0.45 ± 0.04 
for XO-3. VandenBerg & Clem (2003) determined 
empirical color-temperature relations for stars in- 



cluding variations with gravity. We downloaded 
their high temperature table (mentioned in foot- 
note 1 of their paper) and interpolated in the table 
to the effective temperature, surface gravity, and 
metallicity of XO-3 determined above. Doing so 
gives a predicted {B - V)o = 0.440 for XO-3. As 
discussed below, the lightcurve analysis and the 
derived properties of X0-3b favor a larger stellar 
gravity for XO-3. So, we also interpolate the pre- 
dicted colors to a gravity of logg = 4.2 (a little 
more than 3a higher than the spectroscopically 
derived value). Doing so results in a predicted 
{B — V)o — 0.446. These two predicted values 
are within 0.01 and 0.001 mag of the measured 
B — V of XO-3, with an uncertainty of 0.04 in the 
measured color. We therefore find no compelling 
evidence for significant reddening to XO-3 and we 
adopt V — 9.80 ± 0.03 as the extinction corrected 
magnitude of the star. 

We used our spectroscopic effective tempera- 
ture, spectroscopic gravity, and an assumed dis- 
tance to derive a bolometric correction by interpo- 
lating the "high temperature" table from Vanden- 
Berg & Clem (2003). We combined the bolometric 
correction with the observed V-band magnitude 
to determine stellar luminosity. Then we used the 
stellar luminosity and our spectroscopic effective 
temperature, iron abundance, and alpha-element 
enrichment to interpolate the Y^ isochrones to 
produce the probability distribution functions in 
Figure [T] 

The parameters we derive place XO-3 above the 
zero-age main sequence. As a result, reasonable 
fits to the data are possible for pre-main sequence 
evolutionary tracks. However, if XO-3 were indeed 
a pre-main sequence star, we would expect strong 
Li I absorption at 6708 A. For example, in their 
survey of Pleiades F stars, Boesgaard et al. (1988) 
show that the Li I line at 6708 A is approximately 
as strong as the nearby Ca I line at 6718 A is 
stars of similar spectral type to XO-3. We show 
in Figure [8] the Li line region in XO-3 produced by 
averaging the spectra obtained at the HET (this 
region is free of iodine absorption). The position 
of the Li I line is marked and the nearby Ca I line 
is clearly visible. The Li line is not detectable, 
indicating that XO-3 is substantially older than 
the Pleiades, therefore we conclude XO-3 is not a 
pre-main sequence star. 

The most probable distance, mass, radius, and 
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age for XO-3 are 260 pc, 1.41 Mq, 2.13 Rq, and 
2.69 Gyr respectively, for our best fit gravity of 
logg — 3.95. Combined with our measured semi- 
amplitude K — 1471±48 m s"'^, this stellar mass 
implies a mass of MpSini =13.02±0.64 Mj for the 
planet X0-3b. If we assume that XO-3 is at 240 pc 
instead, then the most probable mass, radius, age, 
and gravity are 1.36 Mq, 1.95 Rq, 2.78 Gyr, and 
logg = 4.00, still placing the star well above the 
main sequence. As discussed below, one of the un- 
usual properties of X0-3b is the relatively large ra- 
dius we derive for the planet. The planetary radius 
is approximately proportional to the inferred stel- 
lar radius. The planetary radius could be smaller 
if we could justify a smaller stellar radius; how- 
ever, doing so begins to bring the inferred gravity 
in conflict with the spectroscopically derived value 
oilogg = 3.95±0.062. Taking a 3ct upper limit, the 
largest allowed value for the gravity is logg = 4.14, 
corresponding to a stellar mass of 1.27 Mq and a 
stellar radius of 1.62 Kq for an inferred distance of 
200 pc. Together with the measured velocity semi- 
amplitude, these values imply a planetary mass of 
MpSini = 12.14 ± 0.40 Mj where the uncertainty 
here is just that resulting from the radial velocity 
fit. Clearly, a trigonometric parallax measurement 
for XO-3 will be of great value in better establish- 
ing both the stellar and planetary parameters. 

3.4. Light Curve Modeling and the Plan- 
etary Radius 

For additional photometric analysis and lightcurve 
fitting, all the E.T. photometric data shown in Fig- 
ure [3] were averaged into 7.4 minute bins. Binning 
permitted outlier rejection and empirical estima- 
tion of the noise by the scatter of the individual 
observations, which is helpful in cases such as this 
in which residual calibration errors can be signifi- 
cant. Differences in the depth and the shape of the 
transit light curve in different photometric bands 
are both expected and observed to be substantially 
smaller than the photometric uncertainties in the 
E.T. lightcurves, so we combined the light curve 
from all photometric bands to improve the signal 
to noise ratio. The binning begins with phasing 
the light curves using the refined ephemeris de- 
termined in §3.1. The final binned light curve 
consists of a robust average of E. T. photomet- 
ric measurements in bins 7.4 minutes in duration. 
This robust average is determined by first calcu- 



lating the median and median absolute deviation 
for each bin and rejecting points in the bin that 
deviate by more than 4cr from the median. Next, 
the mean and standard deviation are calculated 
for each bin and points which are more than 4ct 
from the mean are rejected. This rejection based 
on the mean and standard deviation is repeated 
one more time before a final ("robust") mean is 
calculated. The final binned light curve is sup- 
plied in Table [1] and is shown in Figure M along 
with the fits described below. The uncertainty in 
each binned light curve data point is the standard 
deviation of (surviving) measurements in the bin 
divided by the square root of the number of mea- 
surements in that bin. The uncertainties in all 
bins are then multiplied by a correction factor as 
described below (the uncertainties in Table [T] do 
not have this correction factor applied). 

We modeled the mean transit light curve us- 
ing the analytic transit model of Mandel and Agol 
(2002). For a planet with an eccentric orbit, the 
transit model has ten parameters: transit mid- 
point (to), orbital period, stellar mass (M^,), stel- 
lar radius (-R*), planet radius (i?p), inclination (i), 
quadratic limb darkening law coefficients (wi and 
U2), eccentricity (e), and longitude of perihelion 
(w). Throughout the following analysis, the or- 
bital period remains fixed at the period found in 
§3.1. The radial velocity data provide the best 
estimates for e = 0.26 and ui = —15.4 that re- 
main fixed during the minimization. For each 
fit presented below, we fix the stellar parameters 
to specific values based on the spectral synthesis 
and isochrone analysis presented above. We inter- 
polate quadratic limb darkening coefficients from 
Claret (2000) for the Cousins R photometric band- 
pass (ui = 0.220; U2 = 0.380) based on the grav- 
ity, effective temperature, and metallicity ([M/H]) 
determined from the spectroscopic analysis (Ta- 
ble [7]). The remaining parameters, to, Rp, and 
i, are allowed to vary and are solved for using the 
"Amoeba" algorithm mentioned earlier. For input 
into the Mandel and Agol (2002) transit model, we 
determine the projected separation between star 
and planet, S, in the case of nonzero eccentric- 
ity following the procedure outlined by Hilditch 
(2001). For a given i, lj, and e, we determine 
the true anomaly, Ovain, that minimizes 5 (Equa- 
tion 4.9 of Hilditch 2001). The resulting 0min cor- 
responds to the mean anomaly at transit midpoint 
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providing the zeropoint to convert observed times 
to 9 via Kepler's equation. 

We then initially fit the light curve assuming 
the stellar mass (1.41 M0) and radius (2.13 R©) 
derived from the spectral fits as described above. 
The free parameters of the light curve fit are then 
the radius of the planet and the inclination of the 
orbit (and the time of mid-transit). The resulting 
fit is shown in the top panel of Figure [51 We then 
repeated the light curve fits using the extremes 
in the stellar mass (1.33 and 1.90 M0) and ra- 
dius (1.49 and 2.36 R0) which result from the la 
uncertainties in these parameters derived above. 
Together, this then gives us estimates of the plan- 
etary radius (i?p = 1.95 ± 0.16 Rj) and orbital 
inclination (i = 79.°32 ± l.°36) which are also 
given in Table [H The total mass of the planet 
is then 13.25 ±0.64 Mj. Examination of the figure 
shows that the model does not fit all aspects of 
the lightcurve: the model ingress and egress ap- 
pear noticeably longer than the observed ingress 
and egress. 

In order to better model the shape of the ingress 
and egress while still maintaining the total dura- 
tion of the transit, the radius of the star and the 
planet both need to be reduced and the inclination 
must increase so that the total chord length the 
planet traverses accross the face of the star stays 
approximately the same. Figure [7| shows that from 
the isochrone analysis, the stellar mass and radius 
both decrease as the distance is assumed to be 
smaller. We have performed the isochrone analysis 
every 10 pc, so we computed model light curves, 
stepping down in distance 10 pc at a time from 
our favored distance of 260 pc. For each new dis- 
tance, we adopt the corresponding stellar mass 
and radius resulting from the isochrone analysis 
and compute for the fit to the light curve. The 
minimum in occurs between 190 and 200 pc. 
Fitting a parabola to the lowest 4 values, gives 
a best fit to the light curve for a distance of 185 
pc which gives = 1.24 Mq, R.^ = 1.48 Rq, and 
\ogg = 4.19. This fit is shown in the lower panel 
of Figure m As mentioned earlier, the the best 
is quite large (~ 129) for the 63 degrees of free- 
dom in our fit, indicating that the photometric 
uncertainties have been underestimated. We at- 
tempt to correct this by determing the multiplica- 
tive factor required to give a = 63.0 (a reduced 
X^ = 1.0). The uncertainties shown in Figure [9] 



have been multiplied by this factor. Having done 
this, we can then assign as the uncertainty in Rp 
and i the difference between the best fit parame- 
ters and those we obtain when Ax^ = 1.0. Doing 
so, we find that the transit lightcurve is best fit 
for Rp = 1.25 ± 0.15 Rj and i = 83.°32 ± l.°26. 
The total mass of the planet is then 12.03 ± 0.46 
Mj. This value for the planetary radius is con- 
siderably lower than the value derived above from 
the default stellar parameters. As discussed be- 
low, we regard this as a lower limit to the true 
planetary radius; however, we note that the ambi- 
guity described here can be greatly diminished by 
obtaining very accurate photometry of a transit of 
XO-3 and/or by obtaining a precise trigonometric 
parallax measurement. 

4. Discussion 

4.1. Comparison to other Transiting Plan- 
ets 

The ~ 13 Mj planet X0-3b is unusual in many 
respects compared to the sample of known extra- 
solar planets. The mass of X0-3b is quite large 
compared to most of the known extra-solar plan- 
ets. Zucker and Mazeh (2002) noted that the most 
massive short period planets are all found in mul- 
tiple star systems. Udry et al. (2003) emphasized 
the general lack of massive planets on short pe- 
riod orbits, particularly for planets orbiting sin- 
gle stars, and interpreted these results in terms of 
planetary migration scenarios. The star XO-3 is 
not known to be a binary, though it is relatively 
unstudied. While it is in the Tycho-2 catalog (H0g 
et al. 2000), XO-3 does not have a significant 
proper motion measurement from the Hipparcos 
mission. It is not known if any of the nearby stars 
seen in Figure [T] are common proper motion com- 
panions or not. Udry et al. (2003) point out that 
among planets orbiting single stars known at that 
time, there were no planets more massive than 2 
Mj with periods less than 100 d (see also Eggen- 
berger et al. 2004). Using data for 218 extra- 
solar planets as compiled on the Geneva Extra- 
solar Planet website (http://www.exoplanets.eu) 
updated as of 27 May 2007, this general lack of 
short period, very massive planets is still appar- 
ent. There are only 3 other planets (HD 162020b, 
HD 17156b, and HAT-P-2b) orbiting apparently 
single stars with periods less than 30 days and 
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masses larger than 2.5 Mj, and one of those (HD 
162020b) is suspected to actually be a much more 
massive brown dwarf (Udry et al. 2002). There is 
only 1 planet other than X0-3b with a mass larger 
than 2.5 Mj and a period less than 4 days (HD 
120136b, Masim = 4.14 Mj, Butler et al. 1997), 
and this star is known to be in a stellar binary 
system (Hale 1994). 

The eccentricity of X0-3b is also quite rare rela- 
tive to other known extra-solar planets, given the 
short orbital period of X0-3b. The eccentricity 
of the shortest period extra-solar planets are all 
quite low (e < 0.05) and most are consistent with 
zero eccentricity (Halbwachs et al. 2005). This 
is generally believed to be due to either tidal cir- 
cularization (e.g., Wu 2003; Ivanov & Papaloizou 
2007) or the result of smooth orbital migration in 
a disk which is not thought to produce significant 
eccentricity (Murray et al. 1998). The eccentric- 
ity of X0-3b is larger than all other planets with 
periods less than that of X0-3b. We discuss the 
eccentricity of X0-3b further in Section l473l While 
the mass, period, and eccentricity of X0-3b make 
it quite rare among extra-solar planets, it is not to- 
tally unique. The recently discovered HAT-P-2b 
(Bakos et al. 2007) is very similar in many re- 
spects with a period of 5.6 d, Mp = 9.04 Mj, and 
e = 0.52. Also, the hot Neptune GJ 436b has a 
period of 2.6 days and e = 0.15 (Butler et al 2004; 
Deming et al 2007). Apparently the assumption, 
commonly-held prior to 2007, that short-period 
planets should have circular orbits was incorrect. 

It is now well established that host star metal- 
licity correlates with the likelihood of finding a 
Jupiter mass companion orbiting the star in the 
sense that higher metallicity stars are more likely 
to have planets (Santos et al. 2004; Fischer & 
Valenti 2005). The metalhcity of XO-3 ([M/H] = 
-0.20) is relatively low, and the above mentioned 
studies show that planets are found around only 

3% of stars with similar metallicity. There are 
only two other very massive planets {Mp > 5 Mj) 
known around stars (HD 111232, HD 114762) with 
metallicities lower than that of XO-3. The correla- 
tion between host star metallicity and the proba- 
bility of finding a planet is often taken as support 
for the core accretion model (e.g.. Pollack et al. 
1996; Bodenheimer et al. 2000; Hubickyj et al. 
2005) of planet formation. Under that assump- 
tion though, it may be surprising that a few very 



massive planets have also been found around low 
metallicity stars. It has also been suggested that 
host star metallicity decreases on average as the 
planet mass increases which appears inconsistent 
with the core accretion model (Ribas & Miralda- 
Escude 2007). On the other hand, gravitational 
instabilities in a disk are not expected to produce 
a metallicity-planet correlation (Boss 2002). Tak- 
ing this notion further, Ribas and Miralda-Escude 
(2007) have suggested that the sample of extra- 
solar planets consists of objects formed via two dif- 
ferent paths: core accretion in a disk and fragmen- 
tation of a pre-stellar cloud. The high mass, rela- 
tively large eccentricity, and low metallicity of the 
X0-3b system would then "fit" the notion of Ribas 
and Miralda-Escude of a planet formed from the 
collapse of a pre-stellar cloud, but of course the 
observations of X0-3b do not prove (or disprove) 
either formation scenario. 

4.2. How Big is XO-3b? 

Potentially, one of the most unusual and inter- 
esting properties of X0-3b is the large radius of 
the planet. At an age of 2-3 Gyr, giant plan- 
ets and low mass brown dwarfs are expected to 
have radii very close to that of Jupiter (Burrows 
et al. 2001; Baraffe et al. 2003). Using the ta- 
bles in Baraffe et al. (2003), the radius of a 13 
Mj planet is expected to be 1.03 Rj at 1 Gyr, 
and 0.97 Rj at 5 Gyr. Including heating from 
the central star can increase the expected radius 
for so-called "hot Jupiters" (e.g., Bodenheimer et 
al. 2003). Fortney et al. (2007) compute plane- 
tary structure models up to 11.3 Mj for a range of 
orbital separations, including the insolation pro- 
duced by absorption of solar radiation. At an age 
of 300 Myr or older, all the models computed have 
radii less than 1.3 Rj. In all cases analyzed by 
Fortney et al., their models predict a substantially 
smaller radius than the 1.95 ±0.16 Rj determined 
based upon the values of the stellar mass and ra- 
dius inferred from the spectroscopic and isochrone 
analysis of XO-3. Such a large discrepancy be- 
tween the observed and predicted planetary ra- 
dius is reminiscent of TReS-4 which has a radius 
of 1.674±0.094Rj (Mandushev et al. 2007). How- 
ever, the radius of X0-3b is quite uncertain, and 
the light curve itself favors a smaller planetary ra- 
dius of Rp = 1.25 ± 0.15 Rj. Because X0-3b is in 
an eccentric orbit, for comparison with the models 
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of Fortney et al. (2007) we compute the average 
orbital separation between X0-3b and its star and 
then adjust this separation to take into account 
the difference in luminosity between XO-3 and the 
Sun. The smaller planetary radius for X0-3b is 
associated with a smaller stellar radius for XO-3, 
so for this exercise, we take i?* = 1.57 R©. We 
thus compute an effective distance of 0.024 AU to 
use for X0-3b when using the tables of Fortney et 
al. (2007). Interpolating in these tables, the pre- 
dicted radius for X0-3b is 1.18 Rj at 1 Gyr and 
1.11 Rj at 4.5 Gyr. Both values are lower than 
the smaller radius inferred for X0-3b, but only by 
about la. A more precise radius for X0-3b would 
be a very interesting benchmark for the theory of 
extra-solar planets. 

The actual radius of X0-3b can be better con- 
strained by a true parallax measurement and more 
precise photometric observations (ideally in mul- 
tiple colors) of additional transits. The constraint 
from this work that favors the larger radius for 
X0-3b is the stellar gravity inferred from the spec- 
troscopic analysis. A number of investigators have 
pointed out that stellar gravity determinations are 
notoriously difhcult and the relative error in this 
parameter is often substantially larger than almost 
all other measured properties of transiting extra- 
solar planetary systems. Many of these investiga- 
tors have elected to use the transit light curve fits 
in combination with stellar isochrone models to es- 
timate the stellar and planetary radius and hence 
logg for the star (e.g., Sozzetti et al. 2004, 2007; 
O'Donovan et al. 2006; Bakos et al. 2007; Man- 
dushev et al. 2007). The photometry of X0-3b 
is not as precise as that used in those studies, so 
we do not emphasize the photometrically-derived 
parameters as those investigators have. However, 
we can examine the spectroscopically determined 
logg to estimate by how much, and in what sense, 
it is in error. 

The above studies typically use spectroscop- 
ically determined effective temperatures with 
isochrone models to determine the stellar grav- 
ity, implicitly assuming the gravities determined 
from isochrone analyses are more accurate than 
the spectroscopic values. We can compare the 
gravities determined from the isochrones and 
our spectroscopic analysis for stars of similar spec- 
tral type to XO-3. The sample of stars studied by 
VF05 all have accurate Hipparcos parallax mea- 



surements, allowing them to use the isochrones 
with the spectroscopically determined T^ff to de- 
termine the stellar gravity. This gravity can then 
be compared to the gravity VF05 derive from the 
spectra alone, using the same techniques used in 
this paper. Figure [10] shows the result of this 
comparison for 79 stars with Tef / > 6200 K from 
VF05. Recall T^fj = 6429 ± 50 for XO-3. There 
is considerable scatter in this figure; however, in 
the immediate vicinity of the spectroscopic grav- 
ity determined for XO-3, the gravities determined 
from the isochrones are lower. This is the op- 
posite sense to what is suggested by the transit 
light curve of XO-3. As a result, there is no clear 
indication in this sample of stars that the spectro- 
scopically determined gravity for XO-3 is biased 
low. When making this comparison, it is appro- 
priate to consider the accuracy of the gravities 
predicted by the isochrones. Hillenbrand and 
White (2004) report excellent agreement (better 
that 3%) between dynamically determined stellar 
masses from eclipsing binaries and masses deter- 
mined from the Y^ isochrones for main sequence 
stars more massive than ~ 0.6 M©; therefore, we 
find no reason to doubt the Y^ isochrones appro- 
priate for XO-3. We also note that results from 
transit light curve analyses are subject to various 
assumptions and potential biases that may exist 
in the models used for the fit. For example, Auf- 
denberg et al. (2005) confirm the prediction of 
AUende Prieto et al. (2002) that ID atmosphere 
models, including the ATLAS models used by 
Claret (2000) to determine limb darkening coeffi- 
cients, predict too much limb darkening at optical 
wavelengths, depending on the treatment of con- 
vection and convective overshoot in the model. In 
fitting light curves, such a bias will favor models 
with smaller impact parameters and smaller stellar 
radii. Therefore, until better data can be obtained 
for XO-3 (a precise parallax measurement and/or 
more precise transit photometry), the stellar and 
hence planetary radius will remain significantly 
uncertain. At this point, taking the photomet- 
ric lightcurve and the spectroscopic gravity into 
account, we estimate the radius of X0-3b to be 
1.10 < Rp < 2.11 Rj including la uncertainties 
on both limits. The lower limit is almost certainly 
too low however, given the results of Aufdenberg 
et al. (2005) and AUende Prieto et al. (2002) on 
limb darkening showing that the values we have 
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used here are likely overestimated. 

4.3. Tidal Circularization and the Eccen- 
tricity of XO-3b 

As mentioned in Section [O] an interesting as- 
pect of the orbit of X0-3b is its significant eccen- 
tricity, e = 0.260 ± 0.017, because most planets 
with orbital periods similar to X0-3b are believed 
to have already been tidally circularized (Halb- 
wachs et al. 2005). Independent of exactly how 
X0-3b arrived at its current orbit, it is interest- 
ing to consider how long it can remain in such an 
orbit under the influence of tidal circularization. 
The question of tidal circularization of close extra- 
solar planets has been studied by several investi- 
gators (see Adams & Laughlin 2006 and references 
therein) . The equations governing the tidal circu- 
larization of extra-solar planets depend on very 
uncertain planetary quality factors, Qp (for ex- 
ample, see discussion in Gu et al. 2003) and in 
some cases equally uncertain stellar quality fac- 
tors, Q*, if considering tides raised on the star by 
the planet. As an illustration of the difficulty in 
estimating quality factors, Mathieu (1994) point 
out that theoretically determined values of im- 
ply a slow rate for close binary stars to circular- 
ize which is contradicted by observations in stellar 
clusters of various ages. As a result, many in- 
vestigators try to use the observations (of binary 
stars and extra-solar planets) to empirically esti- 
mate the quality factors. For Jupiter mass extra- 
solar planets, Qp - 10^- 10<^ (e.g., Gu et al. 2003; 
Adams & Laughlin 2006). 

Circularization timescales depend linearly on 
Qp, so we assume Qp = 10^ to get predictions 
at the long end of what can currently be esti- 
mated. Using equation (18) of Gu et al. (2003) 
and equation (3) of Adams & Laughhn (2006) 
to estimate the circularization time for X0-3b, 
we find timescales of 0.29 Gyr and 0.33 Gyr, re- 
spectively, using values corresponding to the large 
radius (1.9 Rj) for X0-3b. The circularization 
timescale in both studies depends on the recipro- 
cal of the planetary radius to the fifth power, so 
the circularization timescale grows to 2.92 Gyr and 
3.32 Gyr, respectively, using values corresponding 
to the small radius (1.25 Rj) for X0-3b. Clearly, 
there is considerable uncertainty in the value of the 
circularization time; however, because circumstel- 
lar disks appear to be lost after 10-20 Myr (e.g.. 



Haisch et al. 2001; Mamajek et al. 2004), it ap- 
pears that X0-3b perhaps should have circularized 
by now if it originally arrived at its current loca- 
tion while the circumstellar disk was still in place 
(i.e. through migration). If instead, X0-3b has 
been scattered in to its current location (e.g.. Ford 
& Rasio 2007) more recently, there may not have 
been enough time for circularization to occur, par- 
ticularly for the parameters corresponding to the 
smaller radius for X0-3b. A third possibility is 
that additional, undetected planets orbiting XO-3 
maintain the eccentricity of X0-3b (e.g., Adams 
& Laughlin 2006). It is interesting to note though 
that for planets in short period, eccentric orbits 
such as X0-3b, the tides raised in the planet can 
deposit substantial energy into the planet (e.g., 
Gu et al. 2003; Adams & Laughlin 2006) which 
can inflate it to radii coresponding to the larger 
value we measure (Gu et al. 2003). It therefore 
appears X0-3b can again serve as a very interest- 
ing benchmark for studies of tidal circularization 
and tidal heating once a more accurate radius can 
be established for the planet. 

5. Summary 

X0-3b is a massive planet or low mass brown 
dwarf in a short period (3.2 d), eccentric (e — 0.26) 
orbit, around a somewhat evolved F5 star. There 
is relatively large uncertainty in the planetary pa- 
rameters owing to the unknown distance to XO-3. 
The mass of X0-3b is 11.57 - 13.97 Mj and the 
radius is 1.10 - 2.11 Rj. The larger mass and ra- 
dius are favored by our spectroscopic analysis of 
XO-3, and the smaller values are favored by the 
light curve analysis. A precise trigonometric par- 
allax measurement or more accurate photometric 
light curve date are needed to distinguish between 
these values. 
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Table 1 

XO Survey & E.T. Light Curve Data' 



Heliocentric Julian Date Light Curve Uncertainty Filter Observatory 





[mag] 


(1-cr) [mag] 








2452932.10913 


-0.0017 


0.0027 




1 


XO 


2452932.10938 


-0.0032 


0.0026 




1 


XO 


2452932.11597 


0.0024 


0.0026 




1 


XO 


2452932.11621 


0.0021 


0.0025 


w 


1 


XO 


2452932.12305 


-0.0016 


0.0026 


w 


1 


XO 



''The complete version of this table is in the electronic edition of the Journal. The 
printed edition contains only a sample. 

''Average of N measurements 

'^The filters used in the XO project telescopes are described in McCuUough et al. 
(2005). The bandpass is essentially flat from 4000 - 7000 A. 



Table 2 
All-sky Photometric Magnitudes 



Star'' 


B 


V 


Rc 


Ic 


XO-3 


10.25 


9.80 


9.54 


9.28 


1 


11.37 


10.82 


10.50 


10.18 


2 


12.38 


11.62 


11.17 


10.76 


3 


14.81 


13.07 


12.17 


11.30 


4 


10.13 


9.03 


8.43 


7.88 


5 


12.67 


11.87 


11.43 


10.92 


6 


9.26 


8.78 


8.50 


8.21 


7 


15.52 


13.46 


12.41 


11.50 


8 


15.34 


13.84 


12.59 


11.51 



''Stars are identifled in Figure [T] 
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Table 3 
The Star XO-3 



Parameter 


Value 


Reference 


RA (J2000.0) 


04''21"52^71 


a,b 


Dec (J2000.0) 


+57°49'01".9 


a,b 




9.80±0.03,9.904 ± 0.027,9.86 ± 0.027 


c,b,d 


B-V,Bt- Vt, {B - V)calc 


0.45±0.04,0.451 ± 0.040,0.383 ± 0.034 


c,b,d 


V-Rc 


0.26±0.04 


c 


Rc-Ic 


0.26±0.06, 


c 


J 


9.013 ± 0.029 


e 


J-H 


0.168 d= 0.034 


e 


H-K, 


0.054 ± 0.026 


e 


Spectral Type 


F5V 


c 


d 


260 ± 23 pc 


c 




(-0.6 ± 2.7, 1.6 ± 2.6) mas yr"! 


b 


GSC 


03727-01064 


a 



References. — 

a) SIMBAD 

b) Tycho-2 Catalogue, H0g et al (2000) 

c) this work 

d) Calculated from Tyclio-2 measurements using standard transformations 

e) 2MASS, Skrutskie e al. (2006) 
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Table 4 
Radial Velocity Shifts 



Heliocentric Radial Velocity Uncertainty O — C Bisector 

Julian Date Shift [m s~^] (1 a) [m s~^] [m s^^] Span [m s~^] Telescope 



2454037.0112 


1265 


209 


107 


-31 


HJS 


2454038.0098 


2001 


216 


174 


4 


HJS 


2454137.8215 





134 


-92 





HJS 


2454138.8124 


216 


134 


-112 


170 


HJS 


2454139.8128 


2625 


139 


-117 


133 


HJS 


2454140.8131 


267 


139 


-47 


203 


HJS 


2454141.8047 


77 


149 


13 


112 


HJS 


2454142.7967 


2761 


139 


-1 


68 


HJS 


2454143.7999 


881 


143 


242 


42 


HJS 


2454144.7967 


-35 


161 


34 


152 


HJS 


2454005.8587 


1454 


129 


103 




HET 


2454006.8597 


-976 


122 


16 




HET 


2454113.5842 


-71 


184 


44 




HET 


2454121.7100 


-981 


188 


-43 




HET 


2454121.7190 


-965 


204 


-16 




HET 


2454122.5635 


-1152 


182 


52 




HET 


2454127.7046 


-240 


195 


32 




HET 


2454128.7048 


-1271 


183 


21 




HET 


2454158.6220 


1191 


144 


-32 




HET 


2454159.6091 


-278 


132 


-31 




HET 


2454162.6193 


-30 


194 


-248 




HET 



19 



Table 5 
Transit Timing Measurements 



UT Date 


Observer 


Filter 


Transit Midpoint 
(HJD - 2450000.0) 


2003 Dec 24 


XO 


W" 


2997.72729 


2004 Oct 22 


XO 




3300.90479 


2004 Nov 10 


XO 




3320.05713 


2004 Dec 12 


XO 




3351.99341 


2006 Sep 21 


CF 


R 


3999.86011 


2006 Sep 21 


MF 


R 


3999.85913 


2006 Oct 16 


TV 


R 


4025.39673 


2006 Oct 20 


EM 


V 


4028.58984 


2006 Oct 23 


CF 


B 


4031.77759 


2006 Oct 23 


CF 


R 


4031.77954 


2006 Oct 23 


CF 


V 


4031.78052 


2006 Nov 08 


CF 


B 


4047.73022 


2006 Nov 08 


CF 


R 


4047.73804 


2006 Nov 08 


CF 


V 


4047.73340 


2006 Nov 24 


EM 


I 


4063.69019 


2006 Nov 24 


MF 


I 


4063.69897 


2006 Nov 27 


MF 


I 


4066.88525 


2007 Jan 26 


EM 


I 


4127.53076 


2007 Feb 15 


CF 


V 


4146.67725 


2007 Mar 19 


MF 


R 


4178.58936 



'^The filters used in the XO project telescopes are de- 
scribed in McCuUough et al. (2005). The bandpass is 
essentially flat from 4000 - 7000 A. 
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Table 6 

Orbital Solution and Planetary Parameters X0-3b 



Parameter 


Value 


Notes 


P 


3.1915426±0.00014 d 




tc 


2454025. 3967±0.0038 HJD 




e 


0.260±0.017 




U! 


— 15.4 ± 6.6 degrees 




To 


2454024.7278±0.0570 HJD 




K 


1471±48 m s"^ 




MpSini 


13.02±0.64 


b 


Rp/Rs 


(0.92±0.04) X Rj/Rq 


a 


i 


79.32±1.36 deg 


b,c 


a 


0.0476±0.0005 A.U. 


b 


Mp 


13.25±0.64 Mj 


b,d 


Rp 


1.95±0.16 Rj 


a,b,c 



Note. — 

a) Rj = 71492 km, i.e. the equatorial radius of 
Jupiter 

b) for M, = 1.41±0.08 Mq 

c) for R, = 2.13±0.21 Rq 

d) Mj = 1.8988e27 kg 



Table 7 
Results of the SME Analysis 







Precision 


Parameter 


Mean 


(1 a) 




6429 


50 


logs [cm s ^] 


3.95 


0.062 


V sin i [km s~^] 


18.54 


0.17 


[M/H] 


-0.204 


0.023 


[Na/H] 


-0.346 


0.046 


[Si/H] 


-0.171 


0.017 


[Ti/H] 


-0.262 


0.049 


[Fe/H] 


-0.177 


0.027 


[Ni/H] 


-0.220 


0.033 


[Si/Fe] 


0.006 


0.032 
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Table 8 

Spectroscopically Derived Stellar parameters 



Parameter 


@ 240 pc 


@ 260 pc 


@ 280 pc 




1.33 


1.36 


1.39 


Mass [M0] 


1.36 


1.41 


1.43 




1.39 


1.44 


1.48 




1.92 


2.08 


2.22 


Radius [Rq] 


1.95 


2.13 


2.27 




1.99 


2.17 


2.34 




3.98 


3.93 


3.87 


Log(g) [cm s"^] 


4.00 


3.95 


3.89 




4.02 


3.97 


3.91 




2.66 


2.53 


2.41 


Age [Gyr] 


2.78 


2.69 


2.68 




3.06 


2.83 


2.96 



Note. — For each parameter, the middle row is the 
maximum likelihood value, and the values in the rows 
above and below span the 68% likelihood of the prob- 
ability distributions (cf. Figure [7]). The three columns 
correspond to three assumed distances for XO-3. 
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Fig. 1. — XO-3 is centered, indicated by the two 
hash marks. Stars from Tableware circled. North 
is up; East to the left. The DSS image, digitized 
from a POSSII-F plate with a IIIaF emulsion and 
an RG610 filter, subtends 15' of declination. 




1 .005 



,000 — 



0,995 — 



0,990 — 



0,985 — 



2006.1 1.08 I 
2006.1 1.08 CF R 
"2006.11.08 CF B 
2006.10.23 CF >, 
2006.10.23 CF R 
"2006.10.23 CF B 

_ 2006. 10. 16 TV R 
2006.09.21 MF R 
2006.09.21 CF R 
I I I I Ill 




2007.03.19 MF R 
2007.02.15 CF V 
2007.01.26 EM 
2006.1 1 .27 MF 
2006.1 1.24 MF 
2006.1 1 .24 EM 
I I 



-0.2 



-0. 



0.0 

A t (day) 



0.1 



0.7 



Fig. 3. — Time series photometry of XO-3 from 
2005 - 2007, with dates, observers, and filters in- 
dicated. The observations have been averaged in 
0.006-day bins. The figure is in color in the elec- 
tronic edition. 



Fig. 2. — A total of 2969 individual observations 
of XO-3 by the two XO cameras over two seasons 
2004 and 2005 are shown wrapped and phased ac- 
cording to the transit ephemeris and averaged in 
0.01-day bins (line). The top panel is the full 
lightcurve, and the bottom panel shows the re- 
gion around phase (the primary transit, marked 
with a vertical dashed line in the top panel) en- 
larged. From these data we identified the star as a 
candidate for more refined photometry with other 
telescopes at epochs of expected transits (Figure 
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Fig. 4. — Top: The measured radial velocity data 
are shown along with the best fit orbit. The filled 
circles are based on data from the 2.7 m Harlan J. 
Smith Telescope, and the filled squares come from 
data taken with the 11-m HET. The radial velocity 
curve of XO-3 traces out an eccentric orbit with a 
velocity amplitude K = 1471±48 m s^^, implying 
X0-3b's mass is Mpsmi= 13.02±0.64 Mj. The pe- 
riod and phase used to fold the measured radial ve- 
locities are fixed at values determined by the pho- 
tometric transits. The measured center of mass ve- 
locity with respect to the solar system's barycenter 
has been subtracted (its value is arbitrary since it 
is largely based on the intrinsically relative radial 
velocity measurements from the HJS). Middle: A 
plot of the radial velocity residuals, again phase 
folded with the ephemeris determined from the 
photometric transits. Bottom: The radial velocity 
residuals are shown again, this time as a function 
of time, showing no long term trend. 
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Fig. 5. — The bisector span is plotted as a func- 
tion of the measured radial velocity for all the H JS 
spectra of XO-3. There is no correlation of the 
span with radial velocity and the full amplitude of 
the span variations is less than a tenth of the full 
amplitude of the radial velocity variations. 
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Fig. 6. — The mean spectrum of XO-3 as observed 
(black histogram) and modeled with SME (curve, 
colored blue in the electronic edition) in the re- 
gion of the Mg b triplet. Labels note the elements 
responsible for the indicated spectral lines. Inter- 
mittent line segments (tan) beneath the horizon- 
tal axis indicate wavelength intervals used to con- 
strain the spectroscopic parameters. Very short 
and intermittent line segments (black) immedi- 
ately above the spectrum indicate wavelength in- 
tervals used to constrain the continuum fit. A 
color version of this figure appears in the electronic 
version of the paper. 
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Fig. 7. — The stellar mass, radius, gravity, and age 
which result from the isochrone analysis. The solid 
line shows the derived values as a function of the 
assumed distance, and the gray region shows the 
68% confidence limit on these parameters. Dashed 
lines are drawn at 260 pc corresponding to the 
spectroscopically determined logg = 3.95. 
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Fig. 9. — The combined E.T. transit light curve 
of XO-3 and the best fit transit light curve models 
computed under difi'erent assumptions. The up- 
per panel (a) fixes the stellar mass and radius to 
the values determined from the combined spectro- 
scopic and isochrone analysis. The lower panel (b) 
is determined by minimizing in the light curve 
fit by letting the stellar mass and radius vary ac- 
cording to the isochrone results for different as- 
sumed distances. This best fit produces a stellar 
gravity (logg) more than Scr greater than the mea- 
sured value from the spectroscopic analysis. 
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Fig. 8. — The spectrum of XO-3 in the wavelength 
region of the Li I doublet at 6708 A. The posi- 
tion of the Li feature is marked. The Li feature 
is expected to be approximately as strong as the 
Ca I line at 6718 A in Pleiades age stars of similar 
spectral type, and even stronger at younger ages. 
Therefore, we conclude XO-3 is not a pre-main se- 
quence star, but instead is likely on the post-main 
sequence. 
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Fig. 10. — Comparison of spectroscopically de- 
rived gravities to those derived from isochrone 
analysis for the sample of stars from VF05 which 
have effective temperatures similar to XO-3. Stars 
with [M/H]< 0.04 are shown as blue squares and 
those with [M/H]> 0.04 are shown with red cir- 
cles. The green dashed line in Figure 10 is the line 
of equality. The color version of the figure appears 
in the electronic version of the paper. 
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